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Abstract 



We summarize the signals at HERA in supersymmetric models with explicitly 
broken R-parity. As the most promising case, we consider in detail the resonant 
production of single squarks through an operator L\QiDj, a production process 
analogous to that for leptoquarks. However, the dominant decay of the squark 
to a quark and a photino leads to a very different experimental signature. We 
examine in particular the case where the photino decays to a positron and two 
CN) ■ quarks. Using a detailed Monte-Carlo procedure we obtain a discovery limit in 

the squark mass — Yukawa coupling plane. HERA can discover a squark for a 
| mass as large as 270 GeV and for an R-parity violating Yukawa coupling as small 

ON ■ as 5.8 x 1(T 3 . 

23: 

Oh! 

p; 1 Introduction 

cu : 

Supersymmetry (SUSY) is a promising solution to the gauge-hierarchy problem 
0. It has been thoroughly analysed in the minimal supersymmetric standard model 
^ (MSSM) , which has the minimal field content and the minimal superpotential consistent 

with the Standard Model (SM). The MSSM has a discrete, multiplicative symmetry 
known as R-parity 

R P = (-) 2 *+ 3 ^, (1) 

where S, B, and L are the spin, the baryon and the lepton quantum numbers, respec- 
tively. The fields of the SM, including the extra scalar Higgs doublet, have R p = +1; 
their superpartners have R p = —1. Recently there has been much theoretical 0, cos- 
mological [f|, |5|, |!| [7| and phenomenological work |8], |], |TD, |TTJ on the more general 



supersymmetric extension of the SM with explicitly broken R-parity (R p ). As in the 
MSSM, this only requires the minimal field content consistent with the SM, but in 
contrast allows for all gauge and supersymmetric invariant terms in the superpotential 
including now R p terms which violate either lepton- or baryon-number. R p has been 



shown to be theoretically consistent with grand unification theories [O, fOJ, as well as 



with superstring theory |L4]|, and it is possibly better theoretically motivated than the 



MSSM 0. There has also been an extensive discussion on cosmological bounds on all 
^p-operators §, g 0. In Ref.0 it was shown that most bounds no longer apply, given 
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reasonable assumptions on the mechanism for mattergenesis. We therefore find it very 
much worthwhile to experimentally search for fi p . 

To date the most stringent experimental bounds on the ]/l p Yukawa couplings have 
been determined from virtual effects of the new supersymmetric states in processes 
involving only SM final states. fi p could be found, or these bounds substantially 
improved, by direct searches for supersymmetric particles at colliders. 

The single squark production mechanism is analogous to that of leptoquarks. Thus 
it is important to discriminate fi p squark from leptoquark production. We study in 
detail the distinguishing dominant cascade decay of the squarks to a quark and a photino 
and show that searching for this decay leads to the greatest discovery potential in the 
squark mass- Yukawa coupling plane. 

In the following section we will give a brief theoretical overview of fi, p models. 
Then, in section 3 we will procede to discuss the prospects for discovering supersym- 
metry at HERA, both for the MSSM and for fl p models. We show that, while HERA 
is conventionally considered to be relatively uncompetetive with present limits for su- 



persymmetric searches in the MSSM [[15], [16], [L7|] , it is well suited to study a subclass 
of -^-models, which allow for resonant single squark production. Identifying the most 
promising case, in section 4 we then present cross-section and branching fraction calcu- 
lations. Using these as input to a detailed Monte-Carlo of these processes, and coupling 
this with a simulation of the Zeus detector, we examine possible backgrounds to these 
processes and evaluate the discovery potential in Zeus. We show that fi p leads to 
a much larger mass reach than in MSSM searches and gives HERA the possibility of 
finding supersymmetry. 



2 Theoretical Overview 

Given the particle content of the SM, the most general gauge and supersymmetry 
invariant superpotential is given by [T8 

W = W MSSM + Wf tp , (2) 

where the first term includes those couplings which give mass to the matter fields and 
which are included in the MSSM. The second term in Eq.(^) consists of the J^ p Yukawa 
couplings]] 

W K = X^lLiLjE^p + X' i:jk [LiQjD k ) F + X'^ODj I ),,.-. (3) 

Here we have used the superfield notation; L and E (Q and U, D) are the (left-handed) 
lepton doublet and the antilepton singlet (quark doublet and antiquark singlet) chiral 
superfields respectively. A, A' and A" are dimensionless coupling constants and i, j, k are 
generation indices. SU{2) and SU (3) indices have been suppressed here and throughout. 
The subscript F denotes the supersymmetric invariant part in the product of chiral 
superfields. 

As well as breaking R-Parity, the first two sets of terms in Eq.(|3]) violate lepton 
number conservation, and the last set violates baryon number conservation. Together 

^or a review of supersymmetry and in particular the notation employed here see Rcf. Q. 
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they can lead to an unacceptable level of proton decay, a problem which is avoided if 
R p conservation is imposed (WV = 0). However, it is sufficient to impose a discrete 
symmetry which allows for either only the lepton-number violating terms (LLE, LQD) 
or for only the baryon number violating terms (UDD) in the Lagrangian.^ These 
discrete symmetries are theoretically on an equal footing with each other and with R p 
with the one allowing for the lepton-number violating ^-operators possibly preferred 

!• 

At HERA we are particularly interested in the case where an operator [LiQiDj]p 
is dominant, since this leads to resonant squark production. We therefore present the 
Yukawa couplings explicitly in terms of the 4-component Dirac states 

^-LxQiDj = Kij [-ZLU l L d 3 R - e L u l L d R - {e L ) c u % L d 3 R 

+z> L cf L 4 + v L it L i$ R + (z/ L ) c 44] + h.c (4) 

Here the tilde denotes the scalar superpartners and the superscript c denotes the charge 
conjugated spinor. The minus signs are due to the SU (2) structure of the term. 



3 Super symmetry at HERA 

In discussing the possible R p signals at HERA we shall closely follow the discussion 
and notation of Ref. ||19|| . In particular we shall also make the following simplifying 
assumptions 

• We consider only those cases where in turn one of the 45 operators in Eq.(|3]) 
dominates and the others are negligible, i.e. there is a hierarchy in the Jfi p Yukawa 
couplings. This is a reasonable assumption since the known Yukawa couplings to 
the Higgs in the SM show such a hierarchy (m e i ec /m t0 p — 3 ■ 1CT 6 ). 

• The lightest supersymmetric particle (LSP) is a neutralino, which we denote by 
A. This is a non-trivial assumption, since the cosmological bounds, which restrict 
the LSP to be a neutralino in the MSSM |20j no longer apply in the case of ]fi p 0. 



However, renormalization group analyses in minimal supergravity models, which 
shouldn't be much affected by a small JjL v Yukawa coupling, suggest that the LSP 
is indeed a neutralino f2"T||. When studying specific signals in detail we will for 



simplification assume the LSP to be the photino. 

The phenomenology of $ p differs from that of the MSSM in two main aspects: (1) the 
LSP is no longer stable since it is not protected by a symmetry; it can therefore decay 
within the detector; and (2) it is possible to have single production of supersymmetric 
particles, since the final state is no longer restricted to be i? p -even. We shall make use 
of both of these points when discussing the dominant .^-signal at HERA in section |3.3| . 

(1) The LSP can decay via the process shown in Figure 1, where we have circled the 
R p vertex. It will decay either within or outside of the detector. The condition that 

2 As discussed in Ref.JjJ it is sufficient to exclude the simultaneous presence of the [LQD]p terms 
and the [UDD]p terms. 

3 See Ref.[|| for a brief discussion of relic LSP densities in the case of ft p 
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the LSP decays within the detector translates into a lower bound on the ftp- Yukawa 
coupling [n§ 

A % > I.67 • KT^m^/lOO GeV) 2 (10GeV/M K f /2 } (5) 

where 7 is the Lorentz boost factor in the lab- frame, mj is the mass of the scalar 
fermion which enters via the propagator in the LSP decay (cf. Figure 1). Here we have 
made the approximations of massless final state particles and mj 3> M^, and demanded 

that c~fTL,sp ~ lw, where tlsp is the mean lifetime of the LSP. The constraint (^) is 
always satisfied for the ranges of the squark mass and the R p - Yukawa coupling which 
HERA can probe through single squark production, provided ~ 10 GeV. 

(2) In ]fi p models there are two dominant production mechanisms for supersymmetric 
particles. First, there is the conventional pair or associated production through MSSM 
(super-) gauge couplings. In addition, there is also the possibility of single scalar fermion 
production directly through .^-couplings, which is not possible in the MSSM. In both 
cases the supersymmetric states can decay within the detector to lower mass, i? p -even 
states giving the two sets of signals 

(SUSY Pair production) with {Jft p — decay), (6) 
(ftp, single SUSY production) with (Jjt p — decay), (7) 

which we have summarized in Table 1. The table is analogous to the table presented 



in Ref. [0; for completeness we discuss it in detail in the appendix. Throughout we 
assume that the inequality (|5p is satisfied and thus assume all final states to be R p even. 
If the LSP does not decay within the detector, the signal (Q) reduces to the standard 
scenario in the MSSM with missing px {pr) signals which have been widely discussed. 
For a long-lived LSP the process (^) could still lead to novel signals. However, the cross 
section for single SUSY production is proportional to (A') 2 and thus strongly suppressed 
in the case where the Yukawa coupling A' violates the inequality (|]). In section |3]2| we 
discuss the case (||) and in the remainder of the paper we discuss the more promising 
case (0) in detail. However, before moving on we review the possibilities for observing 
supersymmetry in the case of the MSSM (section |3.1|) . 



3.1 MSSM at HERA 

In the case of conserved R p , the LSP is stable. At HERA, the two dominant production 



mechanisms of the MSSM are ]n| [T6|, [L7J 



e +q -> u e + q, (8) 
e~ + q — > e~ + q. (9) 

In the reaction (^), the decay products are neutral in most scenarios and therefore hard 
to detect. We focus on the neutral current process (0); for an integrated luminosity of 
200 the discovery potential is flCfl PI 



mg + nig ~ 200 GeV. (10) 



4 A more optimistic discovery reach was optained in Rcf . |23|| , however not within the MSSM since 
the neutralino masses were varied independently of each other. 
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In obtaining this limit the authors assumed that there are no decays of the squark or 
selectron other than that to the LSP. 

The discovery reach is to be compared with the existing experimental scalar fermion 
mass bounds. From LEP we have the model independent bounds 



mg, Too ~ 45 GeV, 



which are determined from the Z° width. CDF quote a much stronger, but model 
dependent bound for the squark mass |ES] 



Ttiq ~ 150 GeV, 



(12) 



for which it was also assumed that there are no cascade decays of the squarks and 
that all six squarks are degenerate in mass. In contrast, if one assumes a spectrum of 
non-degenerate squark masses, with one squark substantially lighter than the others, 
the bound on the lightest squark is much weaker [OBJ, nig ~ 64 GeV. Comparing the 



experimental bounds with the discovery reach at HERA ( |10|) we can conclude that for 6 
degenerate squarks there is only a very small window of discovery for the MSSM. If the 
top squark is much lighter than the other squarks |27| , then there is still the possibility 



of discovering R p conserving supersymmetry at HERA |28 . 



3.2 Pair Production with ftp Decay 

In ^models, the LSP decays and the standard j/t signal is either diluted or completely 
absent, depending on the dominant /? p -operator. The decays of the supersymmetric 
particles in the various ^? p -models have been discussed in detail in Ref . 1 19|| . We have 



made use of these results in compiling the dominant signals of the processes (y) in the 
first three rows of Table 1. 

For a dominant operator [LiLjE^F, for example, the LSP decays to two charged 
leptons and a neutrino, 

A -»• ef + el + uj, [L t LjE k } F , (13) 

which in the case of the charged current production (|J) results in a signal of 5 charged 
leptons! (the entry in the first row and the second column). In the case of a dominant 
operator [LiQjD k ] F we have 

A -> (ef, v % , Vi) + 2 jets, [LiQjD k } F , (14) 

and the optimal signal still contains 3 charged leptons. 

From Table 1 we conclude, that in the case of pair production and /^-decays, 
the dominant operator [UDD] F does not lead to any promising signals. The two most 
promising signals for discovery of SUSY are the charged current production (g) followed 
by the LSP decays (0) or (|TJ). 

In order to get a rough estimate of the discovery potential at HERA for these two 
models we make use of the total production cross section a(e~P — > eTq + X) taken 



from reference |L6]. For the dominant operator [LiLjEj^p we restrict ourselves to first 



5 



and second generation charged leptons: i,j, k = 1,2; there is no SM background for 
the production of five charged leptons in the final state. As an optimistic estimate, we 
demand a discovery limit of five events, before any acceptance cuts, for an integrated 
luminosity of 200 pb' 1 to be sufficient for detection. Taking the most optimistic scenario 



for the neutralino sector of Figure 2, Ref [16|, we thus obtain a maximal reach of 



me + mq < 216GeV, [LiLjE k ], i, j, k = 1, 2. (15) 

Next we consider the case (0). For three charged leptons (e or /i), the production rate is 
suppressed by a branching fraction squared (0.88) 2 (cf. Section |4~2j). If we furthermore 
require like-sign tri-leptons the production rate is instead reduced by a factor of (0.44) 2 . 
Taking this latter case, which eliminates all SM background, and again optimistically 
only requiring five observed events for an integrated luminosity of 200 pb~ x we obtain 
as a rough estimate of the discovery limitf] 

mg + nig < 165 GeV, LiQjD k , i ^ 3. (16) 

We now compare these discovery limits with the experimental mass bounds. The model 
independent LEP bounds (|TT|) still hold in the case of fl p . The CDF bounds on the 
squark masses are not valid since they relied on p*T signals. However, using the CDF 
di-lepton search data an independent stringent bound on the squark masses has been 
determined, which incorporates decays of the supersymmetric particles |TD 



m~ q > 175GeV, [LiLjE k ], i,j,k = l,2, (17) 
m~ q > WOGeV, [LiQjDk], i,j ^ 3. (18) 

In deriving these bounds it was assumed that the LSP is the photino and that there are 
five mass degenerate squarks, thereby allowing for the possibility of a light top squark. 
Furthermore it was assumed that the branching fraction BF(q — > q+'y) = 100%. Using 
the CDF di-lepton data the mass bound for a single light squark is as low as the LEP 
bound 0, again assuming BF(q — > q + 7) = 100%. For BF(q — > e ± + q') = 50% and 
BF(q — > q + 7) = 0, which is unlikely, rriq > 108 GeV 0, so that a realistic bound for 
a light squark is presumably somewhat stronger than 45 GeV. 

From our rough estimates we conclude that for five degenerate squarks there is no 
discovery window at HERA for a dominant operator LiLjE^. In the case of an operator 
LiQjDk there is still a small discovery window of about 20 GeV, which warrants a proper 
analysis for i = 2,3 and/or j = 3. For % = 1, the J^ v single squark production at HERA 
leads to a significantly larger discovery potential than the pair production. In sections 
(4) and (5) we discuss this in detail. 

5 We note that for 3 charged leptons in the final state there is the physical background from Z° 
production, e~ + p — > e~ + Z Q + X, followed by the decay Z° — > (ee,fi/i). The total cross section 
multiplied by the relevant branching fraction to two charged leptons (e, fi) is a ~ 2.0 • 10 -2 pb, Ref. p9f . 
About 50% of the scattered electrons vanish down the beam pipe and thus the background is roughly 
a factor five smaller than the signal we are interested in (with no explicit lepton- number violation), 
and should be included in a proper background study. 
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3.3 ^-Production with ^L-Decay 



At HERA, ftp includes the unique possibility of single squark production through the 
couplings given in Eq.(|4]), and the diagram shown in Figure 2a. We have summarized 
the possible signals for single squark production in the last 2 rows of Table 1. Single 
squark production has previously been discussed in Refs. p8|, |30| . In both references only 



the $p-decays of the squarks, as in Figure 2b are studied and in Ref . the discussion 
is limited to the production of top squarks. We extend these analyses to include the 
dominant cascade decays of the squarks, q — ► q + 7, followed by the R p -decay of the 
photino, as in Eq.flTJ]) and shown in Figure 2c. Throughout this discussion we make 
the simplifying assumption that the LSP is the photino, i. e it has no Zino or Higgsino 
admixture. 



3.3.1 Analogies to Leptoquark Production 

The single production of squarks at HERA is analogous to the production of lepto- 
quarks. In the notation of Ref. |3l| the tt, squark corresponds to the R2 leptoquark, and 



the dj squark corresponds to the S\ leptoquark.f] The squarks can decay back to the 
initial state (e~q) through the Yukawa coupling (|4|) just as the leptoquarks do, and as 
shown in Figure 2b. This is observable as a resonant peak in the neutral current deep 
inelastic differential cross section d 2 a(e~p — > e~ X) / dxdQ 2 , Qand we can translate most 



results on the Si and R2 leptoquarks into results on R p . In particular from Ref.pl 
Eq.(8) one expects to be able to observe a down-like squark for m d - ~ 200 GeV and a 
Yukawa coupling as small as A' ~ 0.01. This reach is significantly larger than that of 
Eq. (|16D , where the pair production was followed by LQD mediated decays. 

It is important to discriminate between the discovery of a leptoquark and of super- 
symmetry with ]fi p . There are several suggestions as to how this might be done. 

( 1 ) As can be seen from the couplings in Eq. , it is only the left-handed electrons 
which couple to the squarks, whereas some leptoquarks can couple to both left- and 
right-handed electrons. Thus using polarized beams it could be possible to distinguish 
these two possibilities ptj . However, if the right-handed coupling of the relevant lepto- 



quark vanishes, as it very well may, then this method does not lead to a distinguishing 
signal. 



(2) In Ref.|28| it was suggested that one could distinguish up-like squarks from 
leptoquarks through their different decay modes which are mediated by the Yukawa 
couplings. Through the couplings the U; L can only decay as Ui — > + djR, whereas 
the Si leptoquark has two decay modes Si — > e~ + u, v + d and thus also contributes to 
the charged current deep inelastic scattering. However as we pointed out above, there 
are leptoquarks other than S%. For example, R2 has identical Yukawa-couplings to the 
up-like squark and also decays as R2 — > e~ + djR. Therefore this procedure can not be 
used to distinguish a squark from a leptoquark. 



6 The quantum number F of Table 1 in Ref.[[31j does not apply to the corresponding squarks since 
it was computed assuming the baryon- and lepton-number are conserved at the production vertex. 

7 For the dj case this is also observable in the charged current cross section. 
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(3) In supersymmetric theories the squarks have specific gauge couplings through 
which they can cascade decay to the LSP, q — > q + 7. We study these decays as a 
characteristic feature of squarks and thus study the processes^ 



e + Uj 
e~ + dj 



(djR) 

(UiL) 



djR + 7, 
Ui + 7, 



(19) 
(20) 



which are shown in Figure 2c. The asterisk indicates that the produced squark can be 
virtual. The photino decays as in Eq.(|TJ), and we focus on the decay 7 — > e + + 2 jets. 
This leads to a single positron in the final state, three jets, and no electron or significant 
]/ T . This strongly reduces the background compared to the other photino decays but 



also introduces a branching fraction of about 0.44 (cf. Section |4.2j ). The case of the 
decay 7 — > v e + 2 jets leads to a weak tf T signal with large background, which we do 
not consider promising and have thus also not included in Table 1. 

In the following we assume that the only relevant decays of the produced squarks 
are given by 



d 



el + u iL 
v e + d jR 
djR + 7 
e L + d jR 
u iL + 7 



(21) 
(22) 



which we discuss in more detail in section 14.11. Thus we assume that there are no further 



cascade decays for example to other neutralinos or to the gluino. 



(4) A direct consequence of Eqs.( PT|[2"2"D is that the total decay width of the squarks 
should differ from the corresponding decay widths of the leptoquarks. In neutral current 
deep-inelastic scattering one can in principle measure the mass, coupling and decay 
width of a "leptoquark" resonance. For a true leptoquark the decay width is given in 
terms of the mass and one coupling. If the observed resonance is a squark one should 
find the measured width to be larger than the calculated leptoquark width by up to 
an order of magnitude. However, it is unlikely that the experimental resolution will be 
sufficient to accomplish this p3 . 



3.3.2 Existing Limits 

Before moving on to discuss the processes (|19|,20) in detail we first consider existing 
experimental bounds on the Yukawa couplings A'^-, which lead to resonant squark 
production at HERA. In our study we limit ourselves to the couplings 



1,2; j = 1,2,3. 



(23) 



We do not consider the case where i = 3, as this implies a top quark in the initial and/or 
final state. The case of i = 3 has been discussed in Ref . |p8| , where the authors suggest 



3 To our knowledge this was first proposed in Ref. | 
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that the squarks could be discovered via the leptoquark-analogous ft p -decay. However, 
in that case it is not clear how the squark could be distinguished from a leptoquark. 

The squarks and sleptons can contribute virtually to various processes with i? p -even 
final states through the relevant |? p -operators. The experimental values or bounds of 
these processes directly lead to limits on any contribution beyond that of the SM. In 
our case this leads to bounds on the |? p - Yukawa couplings as a function of the squark 
masses. 

The operators L\Q\Dj, j = 1,2,3, contribute to the semileptonic decays of quarks 



33|| . The (V — A) <g> (V — A) structure of the ftp, effective four-fermion charged current 
interaction is identical to that of the SM. Therefore its contribution is equivalent to a 
shift in the Fermi-constant. If, as we have assumed, one ^-operator is dominant, then 



this is a violation of charged current universality [33]. The experimental limits lead to 



a bound on the corresponding Yukawa coupling at the two sigma level 

^- 003 {m£v)- j = (24) 

and we note that the bound is only correlated to the £77(2) singlet, down-like squark 
mass, not the doublet, up-like squark mass. If we allow for more than one ft v operator 
to contribute, then the violation of universality is reduced and the above bound is 
correspondingly weakened. 

The operators LiQ 2 Dj, j = 1,2,3, contribute to forward-backward asymmetries 
measured in e + e~-collisions. Requiring the combination of ft p and SM contribution 
to agree with the experimental data at the one sigma level one obtains bounds on the 
Yukawa coupling constants |33| 



j * °- 45 Got)' i = 1 ' 2 ' 3 ' (25) 



A' £ 0.26 ( m ^ L ), (26) 
123 V100 GeW v ; 

where the first bound is correlated to the down-like squark mass and the second bound 
to the up-like squark mass. The second bound on A' 123 , Eq. (|26|) , is derived using 6-quark 
asymmetry data and is therefore more stringent. 

The operator L1Q2-D1 also contributes to atomic parity violation, which leads to the 
slightly stricter bound j33| 

depending on the 577(2) doublet up-like squark mass. We see below, that these bounds 
leave a significant discovery window at HERA.[] 

There has also been an estimate of a collider bound on the operators L\Q\D\ which 
could be obtained at CDF [35]. This estimated bound is a function of the slepton 



9 There are further relevant bounds on the couplings A' m from the non-observation of neutrinoless 
double-beta decay Q , and on A' 122 , A' 133 , from bounds on neutrino masses , which are strongly 
model dependent. They are the same magnitude as those given and we omit them in our analysis. 
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mass, not the squark mass and is therefore complimentary to our results. Our bounds 
only very weakly depend on the selectron mass, and we have chosen it to satisfy the 
estimated bound given in Ref. . 



4 Single Squark Production at HERA 

Having identified the most promising process by which ft p might be observed at HERA, 
we now look more closely at this mechanism. 



4.1 Squark Decays 



The squarks can decay to the 2-body final states shown in Eqs.( ]2l"| , |2"2"|) . For the R p 

decays the decay rate is given by 



(A') 2 _ m 2 
167T m 2 



rxK) = ^(i - -3 ) 2 , (28) 



where rriq, m q denote the mass of the initial squark and the final state quark, respec- 
tively. We have neglected the mass of the electron. For the cascade decay to the 
photino, q — > q + 7, we obtain 

/ s e « e2 / m 2 + M? 1/9/ m 2 M? t 
T 2 K = -*-m 9 1 - 9 2 7 /i 1/2 1, -§, -4 29 

where e q is the charge of the squark in units of the electron charge, is the photino 
mass and the phase space function is given by [i(x,y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 
In the limit of massless quarks we obtain for the ratio of the cascade decay to the ftp 
decay 



£2 



> 




(30) 

213, L x QiDj, 

2.8, L x Q 2 D h3 , (31) 
0.95, L X Q 2 D 2 , 

where we have used the Eqs.(E^-E7p to obtain the inequalities. In the case of the 
operators LxQ\Dj we expect the cascade decay to be dominant by more than an order 
of magnitude, whereas in the other cases, due to the weakness of the existing bounds, 
the ftp decay can be comparable in size or even dominant. 
The total decay widths of the squarks are given by 

Ttat(djR) = Ti(m ui ) + T 1 (m dj ) + T 2 (m dj ), (32) 
F to t(u iL ) = T l {m dj ) +T 2 (m ui ). (33) 

The corresponding widths for the leptoquarks Si and R 2 , respectively, are obtained by 
setting T 2 = 0. 
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4.2 Photino Decay 

The photino decays as shown in Figure 1. The differential decay rate for the decay to 
the positron is given by 

c?r(7 -> e + u iL d jR ) aA^ / 2 (7.e)(d.u) 2 (7.11) (e.d) 

dE e dE d ' 2**M^l-my& d V" D l ^ ®l 

2 (7-d)(e.u) _ (7.e)(d.u) - (7-d)(e.u) + (7-u)(e.d) 
+6d Dl CeCu D e D u 

( 7 .d)(e.u) - (7.u)(e.d) + ( 7 .e)(d.u) 

+6d6e D~D d 

(7.d)(e.u) - (7-e)(d.u) + (7-u)(e.d) ^ 

+eu6d DjT d J • (34) 

where a is the electromagnetic coupling, 7, e and d are the 4-momenta of the photino, 
positron and the cij-quark, respectively, and 11 = 7 — e — dis the momentum of the final 
state up-like quark. e e = —1, e u — |, e d = — | are the charges of the final state particles 
in units of the absolute value of the electron charge. The propagators squared are given 
by D e = m?— m\ — 27. e, D d = m| — m|— 27.d + m^, and D u = m? — m\ — 27.u + m^. E e 
and E d are the energy of the positron and the dj quark, respectively. The dot products 
are easily expressed in terms of the masses and E e and E d 

j.e = MyE e , (35) 
7 .d = M^E d , (36) 
7.11 = M|- 7 .(e + d), (37) 

e.d = ~(M? + ml-m 2 d -2j.u), (38) 

and we have neglected the mass of the electron. 

For the other decay modes one obtains analogous formulae. Upon integrating these 
decay rates numerically we obtain the branching fraction for the decay to the positron 

BF(7 - e+ + 2jets) = T (l^ e+ + 2jets) = Q 

"tot (7) 

where we have used = 50 GeV, rriq = 150 GeV. The result only weakly depends on 
these masses. For My = 150 GeV and rriq = 240 GeV we obtain 0.4384. This branching 
fraction was estimated to be \ in Refs. [§, || |10| . If this result (|39D is included in 
Refs. [f|, [l^l the signals are enhanced by a factor of ~ 3.0 leading to stricter bounds for 
example in Eq. (|T7|,[n|) above. 



4.3 Inclusive Reaction Cross Section 

At tree-level the electron-quark differential cross section for the process (|19D is given by 
da(e~Ui — > djty) X' 2 e 2 ( 2 s(s - m 2 ~ - m 2 d ) 2 (i - m 2 d )(i - m 2 ) si 



dtdx 32nP\ 6d Rl +6e Dj + ^ \R S D S \ 
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[u - m 2 d ){u - w 2 ) | | t(fl-m|mg) ' 



where i?| = (s — wi|) 2 + m jr 2 of (cf), Dg = i — rri~, and Z)„ = u — m|. s, i, and 
m are the parton level Mandelstam variables. We have neglected the masses of the 
incoming particles but have included the masses in the final state. The parton level 
cross section for the process given in Eq.(^) can be obtained from the above formula 
through crossing. In order to obtain the electron-proton cross section from the parton 
level cross section we integrate over the product of ( |4~0|) and the appropriate i— th 
generation up-quark structure function Ui(x,T 2 ), 

da(e~P — > dj + 7 + X) r da(e~Ui — > djj) 



dt 



[Z2*pZ».«Ax,T*)dx, (41) 



where x is the usual Bjorken scaling variable and T 2 is the scale of the transfered mo- 
mentum squared. For T 2 we have chosen s, by analogy with photon-gluon fusion (PGF) 
processes. Using instead for instance Mandelstam t for the scale typically changes the 
cross section by less than 5%, though for high squark masses the cross section can be 
enhanced by as much as 15%. In Figure 3 we show the electron-proton cross section 
before integration over x and Mandelstam t (= —Qfjjs), s ° as a function of x at fixed 
t(= 49200 GeV 2 ). 



4.4 Event Rates 

In the MSSM the renormalization group analysis usually leads to five nearly mass de- 
generate squarks, with a possibly lighter sixth top squark. These analyses should not 
be very much affected by a (small) fi p Yukawa coupling. Since we here do not con- 
sider top squark production we set mj. = m„. = M S usy- Furthermore, our results are 
insensitive to to the selectron mass and for simplicity and in order to satisfy the esti- 
mated contraints of Ref . J35| we set mg = Msusy as well. The dominant process is the 
s-channel squark production, so the squark and photino masses will have a significant 
effect on event rates. Figure 4 shows the event rates for a nominal year's running at 
HERA (200 pb^ 1 ) as a function of Msusy for various photino masses. The kinematic 
region considered is; 

A GeV 2 <Q 2 DIS < 9.84- 10 4 GeV 2 , 

1Q- A <x <0.9. (42) 



As determined from Eq. fl4"0|) and seen in figure 3, the parton-level cross section strongly 
peaks at the resonant squark masses and so this is where most of the fi p events are 
produced. The kinematic region (^) includes the fi p squark production peak for prac- 
tically all accessible squark masses and is chosen to minimise theoretical uncertainties 
from DIS physics. The change in gradient in figure 4 at Msusy of around 300 GeV 
is due to the kinematic region used (|4"2] ) - at this point the squark resonant peak falls 
outside the allowed x range. 

Figure 5 shows the cross section as a function of X U j and for fixed rhg, for three 



popular sets of structure function parameterisations |36|, |37|, |38 . 
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One of the main aims of the HERA experiments is to determine parton distributions 
of the proton within the new kinematic regimes avaiable. The present uncertainties 
in the valence quark distributions lie mainly in the low x (< 10~ 2 ) and very high x 
(> 0.4) regions. Using the resonant approximation, m~ pa xs, this corresponds to the 
wide range of squark masses 30 — ► 200 GeV. Therefore in the case of squark production 
from valence quarks (A' Uj -) we find the uncertainties in the production rate due to the 
uncertainties in the structure functions to be small, of the order of a few percent. 
Due to greater uncertainties in the sea quark distributions, uncertainties in the cross 
section for the X' 12 j events can rise to as much as fifty percent, though for most of the 
parameter space they remain at around ten percent or less. These effects are illustrated 
in Figure 5. Plots for varying squark mass were also generated and exhibit a similar 
range of variation for the various parameterisations. In the case of A'^, the cross section 
is heavily dominated by interactions involing the w-valence quarks in the proton, and so 
as can be seen, the differences between the cross sections for j = 1,2,3 are very small. 
The structure function parameterisations were obtained from the PDFLIB package |3T5 



For all subsequent studies in this paper the KMRS BO structure function set will be 
used. 



4.5 Event Topology and Differential Cross Sections 

It is important to study fi p event topologies in order to determine the best cuts to 
separate signal from possible backgrounds at HERA. As discussed in previous sections, 
for the process we are considering 

e + q -> (q')* -> q' + j, ^^e ± + 2jets (43) 

the positron decay channel is seen as the most promising and thus positron transverse 
momentum (pre) is an important quantity in our analysis. Also as we expect both the 
photino and squarks to be relatively massive, thus $ p events generally have large values 
of total transverse energy (Et). Using a Monte-Carlo generator based upon the cross 
section and branching fraction calculations described above, we now look in more detail 
at these events. 

Figure 6 shows plots of da/dpxe and da/dE T for three different photino masses at 
Msusy = 200 GeV. Figure 7 shows the same quantities for different values of M SUSY 
where the photino mass is fixed at 50 GeV. Note the log scale used here, made necessary 
by the strong dependence of the cross section on Msusy- The Jacobian peak at around 
Msusy is evident in the transverse energy plots - particularly for the lower photino 
masses. Also, at lower photino mass a second peak becomes visible at around this 
is the "Jacobian peak" from the photino decay for those events in which the photino is 
emitted with low px- 

It can be seen that pr e and Et are promising quantities to use in selecting SUSY 
events, as was previously shown in Ref. |§. We will also make use of the total energy in 
the detector (Etotal)- The corresponding cuts become potentially more powerful at 
higher SUSY masses, which is also where the cross section falls off and thus efficiency is 
especially important. In the next section we will study some of the possible backgrounds 
to these processes, with particular reference to event identification in the Zeus detector. 
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5 Detection of ]jt p in Zeus 

The Zeus detector is positioned at one of the interaction points of HERA, details of 
which can be found in Ref.^DJ- Here we restrict ourselves to considering the Uranium 



Calorimeter (UCAL) and the Central Tracking Detector (CTD). These are the two 
most important components when searching for high positrons, the most promising 
signature for ]fi p events. 



5.1 Physics Background 

We have considered two types of physics background in detail. The first is heavy quark 
production via photon-gluon fusion (see, for example, Ref . f4T|j ) . These are generally 
Iow-Et events, but they have relatively high cross sections at HERA and thus their 
high Et tail will present a significant background. The most dangerous case can be 
expected to be 6-quark production. The B° thus produced can decay semi-leptonically 
to a positron and a neutrino. The HERWIG generator |^2] was used to generate a sample 
of these events. The cross section for B° production by this process was estimated to 
be 2430 pb with the KMRS BO set of structure functions. This value can vary by as 
much as 30% for other choices of structure function parameterisations, due to the large 
uncertainties in the gluon distribution within the proton. 

The second background considered is W + production via the process 

7* + q -> W + + q, (44) 

where the virtual photon is radiated from the electron and the W + can decay to give a 
positron. Detailed studies of this process at HERA were presented in j|3], |29f . A sample 
of these event was generated with the generator used in the first of these references. The 
cross section for this process multiplied by the branching fraction for the (W + — *> e + u e ) 
was estimated to be 78 fb, varying by as much as 40% for different choices of the 
structure function. 



5.2 'Fake' Background 

'Fake' backgrounds are those caused not by genuine positron production, as in the 
processes above, but by effects inside the Zeus detector faking a high p? positron 
signal. These include 

1. Misidentification of the sign of the curvature (and thus charge) of high pt electrons 
in the tracking detector; 

2. Overlap of hadrons (particularly 7r + ) with photons, providing a positive track in 
the CTD pointing to an electromagnetic cluster in the UCAL. 

In order to investigate the first effect, a sample of high Q 2 deep inelastic scattering 
neutral current events were generated, using the program LEPTO j44|. These events 
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provide a source of high px electrons. The kinematic region over which these events 
were generated was 



300 GeV 2 <Q 2 < 9.84 • 10 4 GeV 2 , 

10" 4 < x < 0.9. (45) 

The cross section within this region is 1350 pb for the KMRS BO set of structure func- 
tions. This value changes by around 10% for the various choices of structure function. 

Studying the second of these effects requires a detailed understanding of the Zeus 
detector and its environment. Previous studies || indicate that these effects can be 
reduced to a level below that of the other backgrounds mentioned. 



5.3 Efficiencies and Resolution 

In order to gain some knowledge of the effectiveness of experimental cuts and thus the 
possible discovery limits for l/t p in Zeus samples of l/t p and background were generated 
corresponding to at least 20 pb" 1 of HERA luminosity. These events were then passed 
through a simple simulation of the Zeus detector, consisting of smearing of particle 
momenta and energies by nominal detector resolutions, followed by geometrical cuts 

17% 



for detector acceptance. The calorimeter resolutions, taken from Ref.KI, wore 



for the electromagnetic calorimeter, and ^7= for hadrons. An angular resolution of 
10 mrad, which is quoted in ||0 as the worst resolution anywhere within the angular 
acceptance of the Zeus calorimeter, was also applied. The CTD resolutions used are 
summarised in Table 2. Where a charged particle entered the CTD, information about 
its angle was taken from the CTD. All energies were taken from the calorimeter. CTD 
momentum information was only used to determine the sign of the charge for positron 
candidates. 

The events were generated at vertices x = y = 0, z = —20cm — > +20cm assuming a 
uniform distribution, where the origin is defined to be the centre of the Zeus detector. 

The ]fi p events were generated for various squark masses and couplings at a photino 
mass of My = 50 GeV, for the photino decays to e~ or e + and two jets. Although we 
attempt to identify only events with positrons in the final state, detector resolution will 
mean that there will be some misidentification of positrons as electrons, and thus some 
mixing between events exhibiting these two types of decay. 

Since the event characteristics change significantly with squark mass, we have de- 
veloped two sets of cuts, one optimised for the lower range of squark masses (we have 
considered squark masses down to 100 GeV) and the other for squark masses at the 
higher limit of observation at HERA. 

For the limit of small X' U j, which is only accessible at low squark mass, the cuts 
rely heavily on the details of the positron signal and must be fairly elaborate. The 
series of cuts found to be optimum in this case, and the corresponding efficiencies for 
Msusy = 100 GeV, My = 50 GeV and A' m = 0.02 are shown in Table 3. Figure 8 
shows the total energy distribution for background and signal before and after these 
cuts, for the same parameters. The cuts which were found to be optimum for the limit 
of high squark mass are summarised in the Table 4. They can be less complex than the 
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cuts for lower squark masses, due to the large amounts of transverse energy present in 
J^ v events at high squark masses. Because the E T cut is very powerful, the cut on the 
transverse momentum of the positron can be relaxed somewhat. Also shown are the 
number of events passing each of these cuts, for M S usy = 240 GeV, My = 50 GeV and 
A' m = 0.07. Figure 9 shows total energy distributions for background and signal events 
before and after these cuts, for M SU sy = 280 GeV, My = 50 GeV and A' m = 0.084 

6 Results 

Figure 10 shows the contour in A'^- space corresponding to five events after 200 pb^ 1 
luminosity and before any detector acceptance or cuts are applied to the data. Also 
shown are the existing limits applying for any or all of these couplings, as detailed in 
sections 4.2 and 5.1. Values to the right of these lines are already ruled out for the 
couplings X'^ as indicated on the plot. 

In order to explain our procedure for obtaining the possible discovery reach at Zeus 
take as an example the most promising case of A' nj -. The existing limit given in Eq . (|24]) 
defines a curve in mass-coupling space. We generated fi p events at various masses 
and couplings spaced along this curve. Applying detector smearing, followed by the 
two sets of cuts from Tables 3 and 4 to these events, we obtain the event numbers 
and signal/background ratios shown in Figure 11. Where the background is so low 
that none of the Monte-Carlo background events passed the cuts, the background value 
corresponding to a one sigma upward fluctuation was used in calculating the ratio. It 
can be clearly seen that the first set of cuts is most effective at the lower squark masses, 
whereas the second set is optimised for masses around 240 GeV. The error bars indicate 
the statistical errors due to the size of the Monte-Carlo samples generated. 

Applying the discovery limit criteria: 

• Five or more Jfi p events surviving our cuts per nominal HERA year (200p6 _1 ). 

• A signal / ^background ratio of at least three. 

to these plots we calculate that should be observable in the Zeus detector up to a 
squark mass of 

Mg ~ 270 GeV, (46) 
for A' iy = 0.081. For M~ q = 100 GeV, % should be observable for 

A' n , ~ 0.0053 (47) 



Carrying out a similar procedure for the 
observable for 



A' 12i cases we calculate that fl p should be 

(48) 
(49) 
(50) 



Mg ~ 205 GeV, A' 121 = 0.53, 
Mg = lOOGeV, A' 121 ~ 0.013, 
Mg ~ 205 GeV, A' 122 = 0.92, 
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Mg = lOOGeV, A' 122 ~ 0.015, (51) 
Mg ~ 175 GeV, A' 123 = 0.46, (52) 



Mg = lOOGeV, A' 123 ~ 0.023. (53) 

It is worth pointing out that it should be possible to extend an analysis by Doncheski 
and Hewett {16] on the effects of virtual leptoquarks at HERA to R p and squarks and 



thus obtain bounds on much higher mass squarks. 

7 Conclusions 

We have presented a specific calculation showing a large expected reach at HERA into 
kinematic regions where R-Parity (ft p ) could be discovered. 

ftp could be discovered for X'nj as low as 5.8 x 10~ 3 or for squark masses as high 
as 270 GeV at X' Uj = 0.081 after the first 200 pb' 1 of HERA luminosity. For the other 
couplings, i.e. A' 12i , the reach in squark mass is somewhat lower, but the weaker existing 
limits on these couplings mean that there is still a substantial window for discovery at 
HERA. These results will hold for photino mass less than squark mass by a few tens of 
GeV or more. 

These limits compare with a possible exclusion limit for Sq leptoquarks after 100 pb~ x 



of 275 GeV at a similar value of coupling A^, obtained in Ref.|45||. We conclude that 
HERA is a very promising experiment for the discovery of supersymmetry with broken 
R p and we eagerly await the results of the present run. 
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A Table of Dominant ^-Signals 

Here we explain in detail Table 1 which summarizes the R p signals at HERA. The 
analogous table for hadron coliders has been given in Ref . [19j and we employ the same 



notation here. We first consider the MSSM pair production mechanisms in the first 
three horizontal boxes, before we come to the single SUSY production via R p in the 
bottom horizontal box. The signals are determined out of the combination of the SUSY 
production mechansim and the ft p decays as indicated in Eqs(|5],|7|). So we in turn discuss 
the production and decay mechanisms. 

At HERA the dominant pair production mechanisms in the MSSM are 

e~ + q — > e~ + q (54) 

e~ + q — ► i> e + q (55) 

as well as 

e - + ( T *,Z°*)^e-+A, (56) 
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where 7*,Z°* denote the virtual neutral gauge bosons. This latter process has been 
discussed in \W7\ for the case of conserved R p . The cross section is relatively low and we 



have not made a numerical estimate for the case of ft v decays. However, for completeness 
we include the dominant signals. 

The respective pair produced SUSY particles e~q, v e q' , and e~A, are listed in the 
first three boxes on the left of the first column. We distinguish the signals according to 
which class the single dominant operator belongs (cf. the assumptions in Section 3). 

The dominant decays of the scalar fermions were determined in Ref . [19| and include 
two different modes. They either decay via a conventional SUSY coupling 



e -> e + A (57) 
q -> q + A, (58) 



or via an ^-operator 



(e, [A, t) + v, LLE 
0(e) = { 2 jets, LQD (59) 

UDD 

LLE 

q -> Q(q) = { (e, v) + 1 jet, LQD (60) 
2 jets UDD 

We have introduced the notation O(f) to collectively denote the decay products of a 
scalar fermion via an ft p operator. The blanks indicate cases, where the sfermion does 
not couple to the operator at tree level. In Ref[0| it was shown that these ft p decays 



are never dominant for operators LLE due to the existing restrictve laboratory bounds 



33|| . For the operators LiQjD^ they can only dominate for higher generational indices. 



The conditions determining which decays are dominant are discussed in detail in |L9 
Clearly the decay products and therefore ultimately the signals depend on the single 
dominant operator under consideration. 

In ftp the LSP also decays within the detector via the diagram shown in Figure 1, 
provided the constraint [| is satisfied. The resulting decay products are 

f e t e l v v LiLjEj, 
A -> <^ (ef , v) + 2 jets, l.,()jl) k .. (61) 
[ 3 jets, UiDjDk, 

where again the decay products depend on the single dominant operator. As discussed 
in Ref. |22|, B5|, in special cases involving only the operators LiQjDj (note the last 



two indices are identical) the LSP can also dominantly decay via the radiative process 
A — » 7 + v. This is discussed in more detail in ||19| . 

The dominant and most promising final state signals for a given dominant operator 
are summarized in the last three columns. In order to determine the complete final 
state we must collect the decay products of all produced SUSY particles. As an example 
consider the first line in the first horizontal box. The production mechanism is e~q — > 
e~q. This is followed by the in most cases dominant cascade decay of the scalar fermions: 

e-q -> e~q + 2A (62) 
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which is the entry in the first column and first row of the first horizontal box. 

For a dominant operator LiLjE^ the LSP decay is given in Eq.floTTj. The index k 
fixes one of the final state charged leptons to be either e, /i, or r. Since the first two 
indices are anti-symmetric, at least one of them is unequal to 3. Thus of the possible 
decay products relating to the first two indices: (e, /x, r) + u, we can choose as the 
promising case: (e, fi) + v. The combined possibilities we write as 

(ee, e/i, pie, \x\x, re, t\x) + v = (e, /i, r) ± <g> (e, /i) + v (63) 

The final state neutrino results in ]/t and since we have two LSPs the total signal is 

e-q+tf r + 2(e,n,T) ± ® (e,/j) T (64) 

which is the entry given in the second column and first row of the first horizontal box 
in the table. 

The other entries in the table are determined by considering the other decay modes 
of the sfermions and the LSP. The second row in the first horizontal box gives the signal 
for the special radiative decay of the LSP: A — > 7 + v. The third row is for the case 
where only the squark decays via a $ p operator and the fourth row of the first horizontal 
box summarizes the signal where both sfermions decay through an lfi p operator. 

The signals for the charged current process e + q — > v + q' as well as the electron- 
photon scattering e + (7*, Z *) — > e + A are summarised in the second and third hori- 
zontal boxes. 

One immediate conclusion is that the operators U DD do not lead to any promising 
signals at HERA. Whereas the operators LLE lead to many charged leptons which 
should make them easy signals to distinguish. 

In the last horizontal box we have summarized the signals for the case of single 
SUSY production at HERA, which is discussed in detail in this paper. The production 
mechanism is the resonant s-channel diagram given in Figure 2a. The squark decays 
are discussed in detail in Section 4.1. The production mechanism is restricted to the 
operators LQD and the final state signals are summarised in column three at the bottom 
of the Table. 



Note Added: Just before completion we obtained a preprint by J. Hewett |48| which 
contains related work. 
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Quantity 


Polar Angle 




90° -> 45° 


45° -> 22° 


22° -> 15° 


1/ Momentum 
(GeV/c)- 1 


0.0018 
+0.0029/p 


0.0012 
+0.0032/p 


0.012 
+0.011/p 


e 


0.0012 


0.00061 


0.0014 


(radians) 


+0.0015/p 


+0.0015/p 


+0.0026/p 





0.00059 


0.00053 


0.0067 


(radians) 


+0.0015/p 


+0.0023/p 


+0.0050/p 



Table 2: Angle and momentum dependence of CTD resolutions used. 



Variable 


Cut 


Events Remaining 


used 




ftp 


bb 


NC DIS 


W + 


None 




1440 


5.1 x 10 5 


2.6 x 10 5 


16 


P Te(GeV/c) 
E T (GeV/c 2 ) 


> 7 

> 40 


537 
527 


185 
139 


123 
123 


12 
10 


E TOT AL(GeV/c 2 ) 


> 40 


527 


139 


113 


8 


MissingE T (GeV/c 2 ) 


< 30 


508 


124 


96 


2 



Table 3: Cuts (1) and Efficiencies, M SUSY = 100 GeV 



Variable 


Cut 


Events Remaining 


used 




ftp 


bb 


NC DIS 


W + 


None 
p T e(GeV/c) 
E T (GeV/c 2 ) 
E T oTAL(GeV/c 2 ) 


> 2 

> 120 

> 450 


152 
63 
56 
56 


5.1 x 10 5 
12700 




2.6 x 10 5 
1550 
26 



16 
13 
0.7 
0.2 



Table 4: Cuts (2) and Efficiencies, M SU sy = 240 GeV 



Figure 1: LSP decay 

Figure 2: a) Single Squark Production b) Leptoquark-like Decay 
c) Cascade Decay 

Figure 3: ft p ep Differential Cross Section 
Figure 4: ]fi p events rates at HERA 
Figure 5: ft p ep Cross Section for different structure function parameterisations 



Figure 6: ft p Differential Cross Sections for various photino masses. 



Figure 7: ft p Differential Cross Sections for various squark masses. 



Figure 8: Total Energy distributions for events containing positrons 
a) Before cuts, and b) after application of cuts in table 3. Msusy = 100 GeV, My = 
50 GeV and A'm = 0.02 



Figure 9: Total Energy distributions for events containing positrons 
a) Before cuts, and b) after application of cuts in table 4. M S usy = 280 GeV, My = 
50 GeV and A' m = 0.084 



Figure 10: ^? p Cross Section contours and existing bounds. 



Figure 11: Signal/Background at discovery limit 
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